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Abstract. This article sketches the proofs of two theorems about sphere packings 
in Euclidean 3-space. The first is K. Bezdek's strong dodecahedral conjecture: the 
surface area of every bounded Voronoi cell in a packing of balls of radius 1 is at 
least that of a regular dodecahedron of inradius 1. The second theorem is L. Fejes 
Toth's contact conjecture, which asserts that in 3-space, any packing of congruent 
balls such that each ball is touched by twelve others consists of hexagonal layers. 
Both proofs are computer assisted. Complete proofs of these theorems appear in 
the author's forthcoming book Dense Sphere Packings. 

1 The strong hexagonal conjecture 

To describe methods, we begin with a proof of the following elementary computer- 
assisted theorem in R^. 

Theorem 1. The perimeter of any bounded Voronoi cell of a packing of congruent balls 
of radius 1 in is at least 4 V3, the perimeter of a regular hexagon with inradius 1. 

If we adopt the convention that the perimeter of an unbounded Voronoi cell is infi- 
nite, then the boundedness hypothesis can be dropped from the theorem. 

Proof. Fix a bounded Voronoi cell of in a packing of congruent balls of radius 1 and fix 
a coordinate system with the center of the Voronoi cell at the origin. The intersection of 
the Voronoi cell with a disk of radius V2 at the origin is a convex disk whose boundary 
C consists of circular arcs and straight line segments. The length of C is no greater than 
the original perimeter of the Voronoi cell. It suflices to show that the length of C is at 
least 4 Vs. 

The boundary C is determined by the set of centers V of balls at distance less than 
Vs from the origin, excluding the ball centered at the origin. 

The following piecewise linear function arises in the proof of the strong dodecahe- 
dral conjecture in three-dimensions. We make repeated use of it. 




(1) 



where /iq - 1.26. 
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Let Ui , U2 e V be distinct points such that T - {0, Ui , U2) has circumradius less than 
^l2. Let ^(ui,U2) be the length of the part of C contained in the convex hull of T, and 
let be the angle at between Ui and U2 (see Figure[T}a).The following inequality has 
been verified by computer using interval arithmetic. 

^(ui , U2) - be{M, , U2) - cL(|ui 1 12) - cL{\M2 1 12)) > 0, (2) 

where b - 4/3 and c - V3/3 - 27r/9 ~ -0.12. Equality holds when T is an equilateral 
triangle with side 2. 

Remark 1. Let Ui(0 and \i2if) be points such that 

|Ul(f)ll = l|U2(0l =2, |Ui(0-U2(0l =f, 

that is a triangle at the origin with sides 2, 2, t. Along this curve with parameter f, the 
constants b and c are the unique choice of constants that give the left-hand side of (|2|i a 
local minimum with value at f = 2. 




Fig. 1. The truncated boundary of the Voronoi cell (shaded) is partitioned into three 
types of pieces, indicated with a thick curve marked I. 



The entire boundary C can be partitioned into finitely many (a) pieces lying in con- 
vex hulls of triangles T, (b) arcs of circles of radius V2, and (c) linear segments from 
u/2, where u 6 V, to a point on the circle of radius V2 centered at the origin (Figure[T]l. 

We extend the inequality (|2]i to boundary arcs of type (b) in the form 

{-bd>0 (3) 

where i is the length of the circular arc and 6 is the subtended angle. This inequality 
is obvious, because { — V2 0, and V2 > b - 4/3. We extend the inequality ^ to 
segments of type (c) in the form 

^(v) - bd{\) - cL{h) > 0, (4) 

where h - || v| /2, is the subtended angle at the origin, and is the length 

of the the segment. When h < ho, inequality (|4]i is a consequence of inequality (|2]l, 
because the segment can be reflected through a mirror to form the two segments in (|2]l. 
When h > ho, the term L{h) is zero. In this case, basic calculus gives the inequality. 
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Let (c be the length of C. We sum the inequalities over the boundary pieces of C of 
types (a), (b), and (c), using inequalities (|3]l, and (|4|i: 

ec-b{2n)-2cY^L{\y\l2)>Q. (5) 

The function L is zero on {/z : h > ho]. We drop such terms from the sum. Lemma[I] 
and Inequality ^ give 

€c > 27rZ7+ 12c = 4V3. 

This proves the theorem. □ 

The following lemma is used in the proof 

Lemma 1. Let V be a set of points contained in a closed annulus at the origin of inner 
radius 2 and outer radius 2ho. Assume that the mutual separation of points in V is at 
least 2. Then 

2l(|v|/2)<6. 

\eV 

Equality is obtained uniquely when V is the set of extreme points of a regular hexagon 
of circumradius 2. 

Proof. In case card(y) < 6, by the inequality L(| v| /2) < 1, it is clear that the sum is at 
most 6, with equality uniquely obtained for the regular hexagon of circumradius 2. An 
easy estimate shows that the angles at the origin between u, v e V is greater than njA, so 
that card(y) < 7. We may therefore assume without loss of generality that card(y) = 7. 

We index the 7 points v,- by / e Z/7Z in their natural cyclic order around the annulus. 
Let 9i be the angle subtended at the origin between v, and v,+i. Let 

Q-/ = arccos(|v||/4) - 7r/6. 

We have the following inequality 

ei>ai + ai+u ieTjlT, (6) 

which is proved by basic calculus: it follows from the intermediate value theorem and 
from an explicit analytic formula for the terms in the inequality I Hall2J . Further, we 
have 

Qf,-0.16L(||v,||/2)-0.32>0 (7) 

which is also proved by basic calculus: by a second derivative test the left-hand side 
is concave as a function of || v, | so that the inequality holds if it holds at the endpoints 
|v, I = 2, 2.52, which is easily checked. 
Summing 0,- over / we get 

27r = 2 > 2 a,. > 2(0.16) ^ L( | V; | /2) + 14(0.32). 
Computing constants, we get 

6>^L(|vJ/2). 

i 

□ 
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Remark 2. The proof can be organized in a way that carries over directly from two 
dimensions to three. In the first step (Lemma[T]i, 2 ^'(0 is shown to be at most the kissing 
number (which is 6 in dimension two and 12 in dimension three). In the second step, the 
estimate of the boundary of the truncated Voronoi cell is transformed into an estimate of 
2 L(-). The second step can be broken into two smaller steps: (a) use a simplex whose 
circumradius is less than V2 to design an inequality with a local minimum at the desired 
solution of to the Voronoi cell problem; (b) extend the inequality from part (a) so that 
it holds on a full geometric partition of the boundary of the truncated Voronoi cell. In a 
final short step, sum all the inequalities to obtain the desired result. 

2 Marchal cells 

In this section, we give details of the partition of the boundary C of the truncated 
Voronoi cell. The partition is based on the partition of Euclidean space into Marchal 
cells l|Mar09l . 




Voronoi cells Rogers simplices k\ composites Marchal A;-cells 



Fig. 2. Partitions of the plane (image source IIHalI2l ). 

Figure|2]shows a packing V of cardinality five. We use the constant V2 to partition 
the plane into levels 0, . . . , 3. Every point has level > 0. For every v G V, we form 
a closed disk of radius V2. A point at level > 1 is a point that lies inside some such 
disk. We form a closed rhombus of side V2 for every pair of distinct points in V whose 
separation is less than VS. By construction, the two points in V are opposite vertices of 
the rhombus. A point of level > 2 is a point that lies some rhombus. We form a closed 
triangle for every triple of distinct points in V whose circumradius is less than V2. A 
point of level 3 is a point that lies inside some such triangle. No point has level > 4. A 
point of level ^ is a point of level > k that does not have level > A: + 1 . 

The points of a given level can be further partitioned using the Rogers partition of 
the plane into simplices |Rog58[ . For each k - Q, . . . ,!>, Rogers simplices that meet the 
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set of level k can be naturally grouped into collections of kl simplices. If P is the union 
of the k\ simplices, then the set of points of level A: in P is called a Marchal k-cell Pk- 

The construction can be generalized to three or more dimensions, again using the 
parameter V2. In n dimensions, the levels extend from to n + 1 in an analogous manner. 
Let S c y be a set of cardinality k+ \ whose circumradius is less than V2. The shapes 
used to define level sets are the convex hulls of 

S UXs 

where Xs is the set of points at equidistance V2 from every point of S . The set Xs 
is sphere of dimension n - card(5). When n — 3, the shapes are balls of radius y/l, 
bicones, bipyramids, and tetrahedra (Figure|3]l. Again in higher dimensions, the Rogers 
simplices can be grouped into collections of k\ simplices, giving Marchal ^-cells P^, at 
each level k. 



card(5) Xs conv(5 U Xs) 



1 


sphere 


ball 


2 


circle 


bi-cone 


3 


pair of points 


bi-pyramid 


4 





simplex 



Fig. 3. Convex hulls used to construct level sets in three dimensions 



Remark 3. Marchal introduced cells to show that the Kepler conjecture in three-dimensions 
can be reduced to an inequality of the form 

2m(|v|)< 12, 

veV 

where M is a certain quartic polynomial, and V is a finite packing contained in a closed 
annulus of inner radius 2 and outer radius VS. 

Remark 4. The book IIHall2l strengthens Marchal's argument to reduce the Kepler con- 
jecture to the inequality 

2l(||v||/2)< 12, (L12) 

\eV 

where L is the function defined above, and V is a packing in the closed annulus with 
inner radius 2 and outer radius 2h(). (In adapting this inequality from dimension two 
to dimension three, the two-dimensional kissing number 6 is replaced with the three- 
dimensional kissing number 12.) The book also gives a computer-assisted proof of the 
inequality ^Lnj , to obtain a new proof of the Kepler conjecture. This article shows how 
to deduce the strong dodecahedral conjecture and Fejes Toth's contact conjecture from 
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Remark 5. An old conjecture by L. Fejes Toth | |Fej72| p. 178] asserts that the minimum 
of 

veV 

is 24 + 14/ V27 X 26.69 as V runs over packings of cardinality 13 contained in a closed 
annulus with inner radius 2 and outer radius VS. The inequality ([T^ gives the best 
known result: 

^ |v| > 24 + 2/!o = 26.52. (9) 

veV 

The inequality also gives upper bounds for the Tammes problem when card(y) = 
13, 14, 15, but these upper bounds are weaker than those known by semi-definite pro- 
gramming LBVOSi . 



3 Strong Dodecahedral Conjecture 

This section sketches a proof of the strong dodecahedral conjecture IBezOOl : 

Theorem 2. The surface area of every bounded Voronoi cell in a packing of balls of 
radius 1 is at least the surface area of a regular dodecahedron of inradius 1. Equality 
is obtained exactly when the bounded Voronoi cell is itself a regular dodecahedron. 



Remark 6. Fejes Toth's classical dodecahedral conjecture | Fej43| is the same state 



ment, replacing surface area with volume. The classical dodecahedral conjecture was 
proved by McLaughlin [HM10 |. To deduce the volume statement from the surface area 
statement, it is enough to use the volume formula Z?/i/3 for a tetrahedron, where B is its 
base area (the face of a Voronoi cell), and /i > 1 is its height. 

Proof. We pick coordinates so that the center of a chosen bounded Voronoi cell is at the 
origin. As in the two-dimensional case, we may intersect the Voronoi cell with a closed 
ball of radius V2. The boundary C after truncation is no greater than before. Let V be 
the set of centers of the packing in the annulus with inner radius 2 and outer radius VS. 

There is a partition C n Pj. of C associated with the set of Marchal A;-cells asso- 
ciated with Rogers simplices at the origin. Write 



area(C) = ^ ai-ea(C n Pk), 



for the areas of the various contributions. Write sol(f t) for the solid angle of the Mar- 
chal cell at the origin, and write dihiPk, v) for the dihedral angle of a Marchal cell Pt 
along the edge through the line through {0, v). 

As a reference cell, we let Po,a be a Marchal 4-cell of the packing giving the regular 
dodecahedron of inradius 1. There exist constants ao and bo > such that 



area(C n + 3flD sol(n) + 3/7D Yj ^(lvl/2)dih(Pi,v) > 0, (10) 
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for all A;-cells Pk and for all V. The constants ao ~ -0.58 1 and bo ~ 0.0232 are uniquely 
determined if we insist that equality is attained when Pk - Pd,a- This inequality has 
been proved by computer by interval arithmetic. 

In more detail, the constants qq and bo are determined by considering a 1 -dimensional 
family PA{t) of tetrahedra with sides 2, 2, 2, t, t, t, for f e R where the three edges of 
length 2 meet at the origin. When t - to ~ 2.1029 (the separation of balls in the ar- 
rangement giving the regular dodecahedron), P^ito) is congruent to Pd,4- Forcing the 
equality to be exact and the derivative to vanish when t = to, we obtain two linear 
equations in two unknowns that determine and bo- 

If we sum ([Tol l over all cells, the solid angles sum to 4n, dihedral angles sum to 2n, 
and L sums to at most 12 by Ineq. ( jLnl l, giving 

area(C) = ^ area(C n Pk) > -3a£,47r - 3i7z)(2;r)(12). 

By the choice of oo and bo, equality is obtained for the boundary Co of the regular 
dodecahedron, 

area(Co) = -3flz547r - 3M27r)(12). 

Hence area(C) > area(C£i). This is the desired conclusion. (The circumradius of the 
regular dodecahedron is less than V2 so that Co is both the truncated and untruncated 
boundary.) Tracing through the case of equality, inequality ( fTOl l is an equality exactly 
when the cell has measure zero or is congruent to Po,4- 



4 Fejes Toth's Contact Conjecture 

L. Fejes Toth conjectured the following result in 1969 | Fej69| , [ Fej89| . 



Theorem 3. In 3-space any packing of equal ball such that each ball is touched by 
twelve others consists of hexagonal layers. 

The proof of this theorem is much longer than the proof of the strong dodecahedral 
conjecture. This section describes the proof strategy. The details are found in lHall2|. 

It is enough to prove that the contact pattern of every ball is the hexagonal-close 
packing (HCP) or face-centered cubic (FCC) kissing arrangement, because these can 
only be extended in hexagonal layers. In fact, the HCP piece has a preferred plane of 
symmetry. Once a single HCP piece occurs, the preferred plane must be filled with 
HCP pieces. A plane forces another hexagonal layer above it and another hexagonal 
layer below it, leading to a packing of hexagonal layers. If no HCP piece occurs, the 
packing is the face-centered cubic packing, which also consists of hexagonal layers. 

Lemma 2. Let V be a packing in which every ball touches twelve others. Then for all 
distinct u, v e V, either |u-v| = 2 or |u-v| > 2/io. 

Proof. Let Ui , . . . , Ui2 be the twelve kissing points around u. Assume that v u,, u. By 
InequaUty ^Ln) , 

12 

L(/!(u, v)) H- 12 = L(/!(u, v)) + L{h{u, Ui)) < 12. 

/=i 

This impUes that L{h(u, v)) < 0, so | u - v | > 2hQ. □ 
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4.1 graph classification problems 

Definition 1. Let be the sphere of radius 2, centered at 0. Let 'V be the set of pack- 
ings V dM? such that 

L cai-d(y) = 12, 

2. VciS^, 

3. ||u-v| e{0,2)U[2.52,4]/orfl«u,vey. 

For each V e 'V, let E^c be the contact graph on vertex set V; that is, the set of 
{u, v) c y such that | u - v | =2. 

Fejes Toth's contact conjecture follows from the Inequality ^Lnj , together with a 
proof that the classification of graphs (y Edc) with y e up to isomorphism contains 
exactly two graphs: the FCC contact graph and the HCP contact graph. 

We formulate Inequality ( |Li2| l as a graph classification problem as well. The in- 
equality holds trivially for a finite packing of cardinality at most 12. For a contradiction, 
we may assume that V belongs to the set of finite packings of cardinality at least 13, 
contained in a closed annulus of radii [2, 2ho] and that violate the inequality: 

2l(|v||/2)> 12. 

Let Esid be the set of edges {u, v) c y such that 2 < |u - v|| < Ih^. The graph classi- 
fication problem equivalent to (L12) is that the set of graphs (y Estd), with V from this 
set of counterexamples, is empty. 

In summary, the proof of Fejes Toth's contact conjecture consists of two graph 
classification problems: one for the contact graphs (y Ecc) involving vertex sets of car- 
dinality 12 and one for the graphs {V,Estd) involving vertex sets of cardinality at least 
13 for the Inequality ( |Li2| l. The proofs of these two classification results differ in detail, 
but the high-level structure is the same in both cases. The graphs are first represented 
combinatorially as hypermaps. (A hypermap is a finite set D together with three permu- 
tations e,n,f of D that satisfy enf = /, the trivial permutation.) A computer program 
classifies the hypermaps satisfying given combinatorial properties obtained from the 
constraints imposed on the graphs (y E). Linear programs eliminate the extraneous 
cases; namely, those hypermaps that exist combinatorially but that do not admit a ge- 
ometric realization. Finally, the inequalities used in the linear programs are proved by 
computer. 

4.2 hypermap classification by computer 

The computer program that classifies hypermaps has been the subject of a exhaustive 
computer code formal verification project by G. Bauer and T. Nipkow IINBS06I . The 
original scope of the project was the set of graphs from the 1998 proof of the Kepler 
conjecture IIHF06I . but in 2010, Nipkow extended this work to include the classification 
of hypermaps needed for for the L12 inequality. 

' There are about 25,000 graphs that arise in the L12 classification and only 8 graphs that arise 
in the contact graph. Because of the vast difference in complexity of these two classification 
problems, our discussion will focus on the L12 classification. 
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4.3 linear programs 

As mentioned above, linear programs eliminate the extraneous cases. The technology 
related to the linear programming has been significantly improved in the years follow- 
ing the proof of the Kepler conjecture. The thesis of S. Obua implemented the formal 
verification of linear programming proof certificates and used this to eliminate about 
92% of the graphs that appear in the original proof of the Kepler conjecture llObuOSI . 
More recent work by Solvyev has optimized the formal verification of linear programs 
to such a degree that the speed of the formal verification of a hnear program rivals 
the speed of the unverified execution of a linear program OSHl II . Work in progress by 
Solovyev intends to make a full formal verification of all linear programs needed to 
prove Inequality ( |Li2| l. 

The linear programs are generated in GLPK from an AMPL model that is indepen- 
dent of the hypermap. An OCaml program generates a separate AMPL data file for each 
linear program. When a single linear program fails to eliminate a hypermap, branch and 
bound methods are used to iteratively subdivide the domain into smaller pieces until lin- 
ear programs are obtained that eliminate the hypermap. The process that was used to 
obtain a system of linear programming inequalities that works uniformly on all hyper- 
maps was fully automated IIHallOI . In brief synopsis, when a linear program fails to 
eliminate a hypermap, two models of corresponding metric graph are compared, one 
based purely on the linear programming estimates of lengths and angles, and a second 
nonlinear model based on nonlinear relations between lengths and angles. A compari- 
son of models is used to determine inadequacies in the linearization. This data is fed to 
a Mathematica program based on various heuristics to construct a candidate nonlinear 
inequality. The inequality is then shipped to the nonlinear optimization package CFSQP 
for extensive nonrigorous testing. From there, a formal specification of the inequality 
is generated in the proof assistant HOL Light. The formal specification is exported to 
program that uses interval arithmetic to verify inequalities by computer; and finally, the 
AMPL model is automatically updated with the new inequality. This process works re- 
markably well in practice to develop a small set of inequaUtie^that uniformly ehminate 
all undesired hypermaps. 

4.4 The classification of contact graphs 

As mentioned in an earlier footnote, there are about 25 thousand graphs that arise in 
the L[2 classification and only 8 graphs that arise in the contact graph classification. For 
the 8 graphs, it was not necessary to follow the lengthy linear programming procedure 
described in the previous subsection. This final subsection sketches a much simpler 
procedure to eliminate the unwanted cases. (Two of the eight possibilities are the HCP 
and FCC, and the other six cases are unwanted.) 

Five of the six are eliminated with linear programming inequalities. The linear pro- 
grams are based on the following simple inequalities: 

1 . The angles around each node sum to In. 

2. The angle of each triangle in the contact graph equals arccos(l/3) ~ 1 .23096. 



About 500 inequalities occur. 
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3. The opposite angles of each rhombus are equal. 

4. Each angle of every rhombus is between 1 .6292 and 2.16672. 

The final case is the graph shown in Figure |4] It is eliminated with the following 
observations. The perimeter of a spherical hexagon with sides 7r/3 is 2n. However, the 
hexagons in the graph are spherically convex, and 27r is a strict upper bound on the 
perimeter of a spherically convex hexagon. Thus, this case does not admit a geometric 
reahzation as a contact graph. Fejes Toth's contact conjecture ensues. □ 




Fig. 4. This planar graph is not a contact graph. 
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